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ABSTRACT 



We compare measurements of horizontal flows on the surface of the Sun using 
hehoseismic time-distance inversions and coherent structure tracking of solar 
granules. Tracking provides 2D horizontal flows on the solar surface, whereas 
the time-distance inversions estimate the full 3-D velocity flows in the shallow 
near-surface layers. Both techniques use HMI observations as an input. We 
flnd good correlations between the various measurements resulting from the two 
techniques. Further, we flnd a good agreement between these measurements and 
the time-averaged Doppler line-of-sight velocity, and also perform sanity checks 
on the vertical flow that resulted from the 3-D time-distance inversion. 



Subject headings: Sun: helioseismology - Sun: interior 
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1. Motivation 

The Sun is a very dynamic system, where plasma flows play a very important role in 
solar dynamo processes. The solar dynamo changes the solar magnetic flelds and gives rise 
to all phenomena of solar activity. Thus the measurement of plasma flows is an important 
part in understanding the physics of the Sun. Plasma flows in the solar photosphere were 
studied by many authors by many methods in the past, some methods are straightforward 
and intuitive whereas some are more complex. Although there is general agreement on 
the multiscale structure of these flows, details still differ with the method used to infer 
these flows. It is not easy to answer the question, which method gives the most accurate 
results. However, by directly comparing various methods one can gain confidence in the 
trustworthiness of the results. 



1.1. Comparisons between methods 

Not counting the direct measurement of the line-of-sight component of the three- 
dimensional plasma flow vector, there are two principal methods for measuring velocity flelds 
on the solar surface. The flrst is based on following structures (granules, supergranules, 
magnetic elements, etc.) carried by the underlying flow in the series of frames capturing the 
same region on the Sun. These provide in general two (horizontal) components of the flow. 
The second (helioseismology) is based on the analysis of the propagation of solar waves 
through the convective envelope and can, in principle, provide us with all three components 
of the flow velocity vector. 



Both groups of methods hav e been carefully teste d, usually by utilising known inputs 



from numerical simulations (e.g. 



Rieutord et al 



the methods to data from different sources. 



2nnih or 



Svanda et al. 



by com paring applications of 



(120061 ) validated and calibrated 
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the Local Correlation Tracking (LCT) method ( iNovemberl 119861) applie d to M ichelson 



Roudier et al 



torn applied 



Doppler Imager (MDI) full-disc Dopplergrams. Recently 
their Coherent Structure Tracking (CST) code to both full-disc Helioseismic and Magnetic 
Imager (HMI) and Hinode intensity grams. The agreement between the results inferred from 
completely different sources was great, except for a systematic error in lower- resolution 
SDO/HMI velocity maps near the solar limb. A correction curve was derived based on this 
comparison. 



Helioseismic meth ods were validated in a similar manner. iGeorgobiani et al.l ( 120071 ) 



applied time-distance (iDuvall et al. 



19931) methods to a numerical simulation of near-surface 
convection, where the solar waves are naturally excited. They used only the surface gravity 
(/) mode to derive the horizontal surface flow field by means of scaling the travel-time 
maps (thus no inversion was performed; such approach is acceptable for the /-mode when 
the surface flow is in question). The comparison of the horizontal [v^ and Vy components) 
flow maps to the known flows in the simulation ended up with correlation coefficients highly 
positive (0.70 for Vx and 0.73 for Vy). At the same time they also applied LCT to the series of 
simulated intensitygrams with a correlation to the known flows of 0.99 for both components. 
However, the authors did not m ention a direct comparison of time-dist ance and LCT 



veloci ty maps. 



Zhao et al. 



Georgobiani et al. 



( 120071 ) used the same numerical simulation as 

(|2007l ) to validate a proper time-distance inversion for flows with a correlation of the 

resulting surface flow maps with the known flows in the simulation of 0.72 for v^, 0.64 for 

Vy, and —0.72 for vertical Vz- Note especially the significant anti-correlation in the vertical 

flow inversion. It was discussed that this is due to the presence of a cross-talk (i.e. leakage 

of the horizontal components into the vertical one in the mass- conserving flow). 

I 



This conclusion was largely confirmed by iSvanda et al.l (1201 ll ) , who validated a 
brand-new time-distance inversion code against synthetic data coming from a numerical 
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simulation of Sun-like convection. Only the inverse modelling part of the whole processing 
pipeline was validated, which allowed avoiding possible problems in the travel-time 
measurements and to study all components of the inversion in detail. It turned out 
that, indeed, for weak vertical flow inversions the cross-talk contributions from horizontal 
components are larger than the proper signal of the vertical flow. The cr oss-talk coraponen t 



was also highly anti-correlated, which explains the negative sign noted by iZhao et al. 



fl2007h . 



A method of minimisation of such cross-talks was invented and successfully implemented. 
The flow estimates from the inversion compared well to the input flows from the numerical 
simulation with a correlation coefficient of 0.9 for horizontal components in the near-surface 
layers and of 0.8 for the vertical component (before the minimisation of the cross-talk it 
was —0.63). 

Recently, similar end-to-end validation of the inversion pipeline for helioseismic 
holography (including t he measurement of helioseismic observables) was performed by 



Dombroski et al. 



( 120131 ). The validation was successful in terms of horizontal components 



and failed (again, due to the cross-talk contributions) for the vertical component. 

Comparisons between different methods for retrieving velocity fields were also done in 

r — 1 1 — I 

the past. iHindman et al.l (|2004)) compared the results of two distinct helioseismic methods, 
ring-diagram and time-distance inversions, with high a correlation (~0.8) between the 
different fiow estimates. A general conc l usioii was that the two distinct helioseismic methods 



observe the same Sun. 



De Rosa et al. 



fl2000l ) compared the horizontal divergence of the 
fiow computed from the horizontal velocity maps measured by means of LCT and a proxy 
for fiow divergence from time-distance /-mode travel-time maps. The correlation of the 
two methods was 0.89 when the different resolutions of the ho rizontal d i verge nce estimates 



was taken into account. Inconclusive results were obtained by lAmbrozl ((20051), using flow 



maps derived from a series of low-resolution full-disc magnetograms measured, by means of 
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LCT when compared to down-sampled velocity maps from time-distance helioseismology. 
The correlation coefficient describing the match of the velocity maps resulted from the 
two methods was close to zero, however, there were compact and continuous regions of 
characteristic size from 30 to 60 heliographic degrees with a good agreement between the 
two methods, so th at one could not conclude that the results were completely different. 



Svanda et al. 



(120071 ) directly compared surface flow maps obtained using LCT applied 
to SOHO/MDI full-disc Dopplergrams and using time-distance helioseismology. A 
map-to-map correlation coefficient was 0.82 for the component of the flow in the direction of 
solar rotation (which was not removed and thus the correlation was dominated by rotation) 
and 0.58 for component in the south-north direction. The effective resolution of maps 
compared was set by the size of the LCT apodization window and was around 44 Mm. 

To our knowledge, the direct validation of the time-distance inversion pipeline 
(including the travel-time measurements) and tracking measurements on smaller-than- 
supergranular scales has not been published. However, it is an important task and must 
be accomplished when the robustness and trustworthiness of helioseismic flow maps comes 
in question. The aim of this study is to perform end-to-end testing of the time-distance 
inversion pipeline running at Max-Planck-Institut fiir Sonnensystemforschung (MPS) in 
Katlenburg-Lindau, Germany and newly implemented also at Astronomical Institute of 
Academy of Sciences in Ondfejov, Czech Republic. By testing the measurements of surface 
flow we inherently shed more confidence on our knowledge of the flow below the surface and 
measured by the time-distance helioseismology, which certainly has consequences to the 
physics of these flows. 
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Measured surface flows 



We measured surface flow fields using the two distinct methods (tracking and time- 
distan ce inversion) applied to data series obse rved by Helioseismic and Magnetic Imager 



(HMI; 



Scherrer et al. 



201 



Schou et al 



2OI2I ) aboard the Solar Dynamics Observatory 



satellite on the 12th and 13th of May, 2010. The results were compared to test the 
correctness of the helioseismic inverse modelling. 



2.1. Flow estimate from granulation tracking 



O ne implementation of tracking methods has a form of coherent str ucture trackinc 



Rieutord et al. 



20071 ). The new implementation of the algorithm (IRoudier et al. 



2OI2I ) splits into three main steps: 



1. In each frame of the series of intensitygrams, individual granules are segmented out. 
The segmentation is based on localisation of local maxima of intensity through the 
intensity curvature, which are tagged as granular centres. The edges or granular 
cells are detected by extension of regions around the detected centres with points 
whose minimal curvature value (evaluated in eight directions from surrounding pixels) 
is greater than a given threshold, while keeping the minimal distance of one pixel 
between each pair of granules. In this way, there is a control over the size of segmented 
structures. Each granule is then identified and tagged. 

2. For each granule, a trajectory throughout the data series is drawn by cross-identifying 
the granules in consecutive frames of data series. Strong assumptions apply in this 
step. Granules are cross-identified taking into account a common granule area in 
the consecutive frames of the series with a threshold that is determined by an upper 
bound on velocity (measured as displacement divided by time lag between the frames). 
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Typically, velocities larger than 7 kms~^ are rejected. The detected trajectories are 
usually averaged over a time window, the average velocity is then computed as total 
displacement change divided by the length of the time window. 

3. It is clear that the velocity field derived in the previous step is irregularly spaced. 
For further analysis, it is necessary to approximate the velocity field by the best 



Mallat 



19891) 



differentiable field. That is done using multi-resolution analysis (e.g. 
where the irregular velocity field is expressed in wavelet components with different 
scales. The differentiable field is then reconstructed from these components. 



The CST algor i thm 



Rieutord et al. 



2008 



las been utilised in many papers (e.g. iTkaczuk et al.l 12007 



2010 



Roudier et al. 



2012 



to name a few) and has proved its usefulness 



in studies of surface velocity fields. 



We applied the code to a series of HMI full-disc intensitygrams. The implementation 
of the CST algorithm used in this study is not sensitive to signal of oscillations in intensity 
images, thus it is not filtered out. The total data series spanned two full days and this was 
split into shorter series to capture the evolution of surface fiows. Each map was computed 
with a 30 minute time window, thus we computed 96 maps of the horizontal components 
of the fiow over the studied two days. The results computed using the full-resolution of 
HMI were binned seven times to incr ease signal-to- noise ratio. The random-error estimate 



is 250 m/s for each 30-minutes map (iRoudier et al. 



20121 ). In these maps, motions within 



supergranules visually dominate the fiow field. 



2.2. Flow estimate from time— distance inversion 



The saying "someone's noise is someo ne else's signal" i s especially true for 
helioseismologists (for a recent review, see iGizon et al.ll2010[ ). The time series of 
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observations contains not only the convective motions but also solar oscillations. In the 
previous methods, the signal of solar oscillations did not carry any useful information that 
could be used by the CST algorithm, and could thus be easily removed without having 
effect on the analysis. However, the solar oscillations observed in surface measurements are 
manifestations of seismic waves travelling through the solar interior, where their propagation 
is affected by anomalies in plasma state parameters. The deviations in travel times of solar 
waves can be measured and used in order to dig out the information about these anomalies. 
Here we focus on the inference of flows, which have the largest impact on difference travel 
times, i.e. the difference between the measured travel times of waves travelling in opposite 
directions. 



Time-distance helioseismology ( iDuvall et al. 



19931 ) is a set of tools useful for measuring 



and interpreting wave travel times. It consists of the following steps: 



1. Wave travel-times are at best measured from a series of full-disc maps of line-of-sight 
velocity (Dopplergrams). HMI provides such measurements with high cadence (45 s) 
thus providing us with almost ideal datasets. The datasets are tracked and remapped 
to Postel projection using standard data processing techniques. We track only the 
disc-centre region (512x512 pixels) with a pixel size of 1.46 Mm always for 12 hours. 
Thus we obtained four consecutive 12-hour Dopplergram datacubes suitable for 
travel-time measurements (see Tab. [1]). We focus on a small patch near the disc 
centre, thus Postel projection approximates the Cartesian coordinate system {x,y,z), 
where x is in the direction of solar rotation, y in direction solar south to north 
and z indicates height. In the small field-of-view, the deviations of Postel-projected 
coordinates from the Cartesian system are negligible. 

2. The datacubes were spatio-temporarily filtered to separate different modes of solar 
waves. In this study we use only the surface gravity (/) mode. 
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3. The travel-times were measured following the approach of iGizon fc Birch! ( 120041 ) . It 
provides linearised fitting of travel-times from measured signal cross-covariances and 
is robust in the presence of noise. Travel-times were m easured for centre-to-annulus 



and centre-to-quadrant geometries (IDuvall et al. 



19971 ) with radii of the annuli 5 to 



20 pixels, thus providing a set of 48 travel-time maps for each of four Dopplergram 
datacubes. 



4. These travel-times are inverted for flows using inversion weights, which is described 
separately. 



2.2.1. Time-distance inversion 



We assume that there is a linear relationship between the measured travel-time 
deviations 6t and velocity vector v given by 

5r'*(r)= [ K'^ir' -r,z)-v{r',z)d^r'dz + n''{r) . (] 

Jq 

Here K"' = { Kt, K^, K'j) is a vector travel-time sensitivity kernel computed in the Born 



approximation (iBirch &: Gizon 



20071 ). r = {x,y) is a horizontal position vector, and index a 



uniquely refers to a particular geometry of the travel-time measurement. 

We attempt to invert for v from ([1]) when knowing 5t and K. A realisation of the 
random noise n"- is not known, however its covaria nce matrix was meas ured from the 
data using the 1/T fitting approach (for details see 



Gizon fc Birch 



2004J ). where T is the 



length of observation. The i nversion can be solved by m eans of a Subtractive Optimally 



Localised Averaging (SOLA; 



Pijpers fc Thompson 



19921) aDProa.ch. T he SOLA technique is 



Jackiewicz et al.ii2008i) . It aims to construct 



standard in time-distance helioseismology (e.g. 
a spatially bound averaging kernel by linearly combining the set of sensitivity kernels while 
keeping the error magnification under control. The estimate for the flow component f J; 



.inv 
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where a = x^y^z, is then given by 

<^(ro, ^o) = J2 [ <(^' - ^o)r'^(r')dV', (2) 
where ly" are the inversion weights to be determined. 

The inverted flow velocity is then a combination of (1) a true velocity component 
smoothed by the averaging kernel /C^, (2) a crosstalk from other components, which we 
attempt to minimise, and (3) a random-noise component f^™*^®, root-mean-square value of 
which {a a) we bound: 

vT{r,,z,)= I lCl{r-r,,z-z^)v^{r,z)d\dz+Y, [ JC^^ir-ro, z; zo)vfs{r, z)d^rdz+vr'%ro, zo] 

(3) 



In terms of the weights, the component of the averaging kernel /C^ is expressed by 

/C^(r, z- zo) = Y. f <(^'; ^o)i^^"(r - r', z)dV'. (4) 

a 

The computation of the inversion weights is a mathematical problem, which usually ends 
up in a ne ed to invert a large ill- posed and ill-conditioned matrix. We use the multichannel 



approach (jjackiewicz et al. 



2012) assuming the spatial invariance of the background. The 



multichannel formulation reduces the computational demands significantly. 

We used the in version code that was already validated using synthetic data 



(ISvanda et al 



20111 ). The inversion we performed was one of the simplest ones. In 
accordance with the travel-time measurements described above, we utilised sensitivity 
kernels for the / mode, where a refers to one of the 48 possibilities combining one 
of three centre-to-annulus and centre-to-quadrant geometries and one of sixteen radii of 
the annuli. The exception was only the inversion for the vertical flow, where we used 
only the centre-to-annulus geometry with the combination of 16 distances. The other two 
geometries, sensitive to waves travelling in the east-west and south-north directions do not 
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contain much information about the vertical flow, do not contribute to the averaging kernels 
and just increase the level of noise. We set a requirement on the averaging kernel so that 
its component in the direction of inversion /C° is strongly localised near the surface (see 
Fig. [1]) with a horizontal FuU-Width-at-Half-Maximum FWHMh = 10 Mm a nd cross-talk 



components /Co for /3 7^ a are minimised (for details how to do so we refer to 



Svanda et al. 



20111). The bounds on root-mean-square value of the random- noise component v'^°^^^ are 
set to be 25 ms~^ and 4 ms~^ for horizontal and vertical components respectively when 
assuming the averaging over 12 hours. 



2.3. Alignment of flow maps from the two methods 

Unfortunately, although the data entering the two analyses originate from the same 
telescope and are thus perfectly co-aligned in the beginning, in the course of processing, the 
coordinate system was changed and the resulting flow maps had to be carefully co-aligned 
again. 

The 2-D flow maps from CST were re-projected onto coordinates grid of the helioseismic 
datacubes utilising Postel projection with pixel size of 1.46 Mm, and averaged in time 
to match time-distance time span. For each of four datacubes, 24 CST flow maps were 
averaged to form one 12-hour average. Averaging over 12 hours increases the signal-to-noise 
ratio again. Assuming that the realisation noise decreases as y/T, where T is the length of 
the time window, the estimate of the noise level in the 12-hour velocities is 50 ms~^. 

Exact coordinates of the projection centres and time span are given in Table [H The 
angular differential rotation is of the form 



u{b) = 2.6373 - 0.3441 sin^ b - 0.5037 sin^ b /irad/s. 



(5) 
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where b stands for the hehographic latitude, was subtracted from re-projected Vx component. 
Using this rate, the datacubes for time-distance analysis were tracked, thus it is natural to 
remove precisely this rotation profile from the non-tracked CST flow estimate. 



3. Comparison of flow estimates 

The co-aligned flow maps may be compared directly on a pixel-to-pixel basis. Both 
methods aim to measure plasma motions on the surface and should ideally match within 
the noise levels. That actually is almost true. Results coming from both methods 
have a different effective resolution. The CST flow estimates were computed using the 
full-resolution HMI intensitygrams and then binned seven times, thus the effective resolution 
is 3.5". The flow estimates from the time-distance inversion are obtained with the averaging 
kernel with FWHMh = 10 Mm, thus having an effective resolution of 13.7", much coarser. 
Thus, the flow maps from CST pipeline must naturally contain flows on smaller scales, 
which is clearly seen e.g. in Fig. [2] when comparing the left and right columns. To make 
flow maps from both methods comparable we convolved the velocity fleld from the CST 
code with a cut through the averaging kernel /C° taken at 2; = 0. The smoothed components 
of the CST flow are also displayed in Fig. [2] in the middle column. 

Pixel-to-pixel statistics for all four datacubes are summarised in Table [21 It gives 
the correlation coefficient CC for both horizontal components (for a curious reader, the 
correlation coefficient with non-smoothed flow estimate from CST is given in parentheses). 

The comparison by eye from Fig. [2] reveals similarities but also some differences. To 
make the comparison more straightforward, we also plot values coming from both methods 



against each other in Fig. [31 The slope of the best 



using random errors in both variables ([York et al. 



east- squares fit through the data points 



20041 ) is also given in Table [21 The 
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values measured by seismology are systematically smaller (by 3-5%) than those measured 
by granule tracking. Such discrepancy can easily be explained by additional averaging 
of helioseismic flow estimates in the vertical direction. When assuming that the flow 
magnitudes peaks just at the surface and its magnitude drops rapidly in larger heights and 
decreases also in depth, the smoothing with the averaging kernels leads naturally to lower 



amplitudes. Possible 
1.6 Mm sugges ted by 
pointed out by 



arge-magnitude (almost 1 kms ^) horizontal flows at depths around 



Duvall &: Hanasogd (|2012| ) do not affect this conclusion, because as 



Svandal (|2012| ). in inversion regularised strongly about the noise term, which 
is also the case here, the possible large-amplitude flows are largely smeared by the extended 
side-lobes of the averaging kernel. The spread of points seen in Fig. |3] seems to be very 
consistent with the levels of random errors estimated for both methods. 

Additional source of differences originates from the fact that both methods may 
not sample (or sample differently) the horizontal velocities of the same layers: granules 
move because their convective flux moves so the sampling function (vertical extend of the 
averaging kernel) of granules is the shape of their convective flux. On the other hand /-mode 
has another vertical proflle which leads to another vertical sampling of the horizontal flow. 

The location of the points, where the difference of the flow components measured by the 
two methods is signiflcant (i.e., the speed values are above their respective noise levels for 
both methods and the difference is higher than the threshold discussed further), is random 
in the fleld of view - these points do not correlate with the respective flow estimates. The 
threshold was determined from a simple model of the expected variance of the difference of 
two random variables, i.e., o"difrerence = \/ ^'Lt + ^Id- The value of the theoretical threshold 
(55 ms~^) is close to the variance of 62 ms~^ determined from studying the histogram of 
the differences of flow estimates from the two methods. The use of the latter value does not 
change the previous statement. 
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The comparison presented here serves as indication, that the hehoseismic pipehnes 
provide reasonable estimates of flows. It turns out that by directly comparing the velocity 
fields coming from granule tracking on the surface, we do not see any significant systematic 
errors. 



4. Vertical flow 

The helioseismic inversion allows us to measure the relative weak vertical flow at 



the surface (thanks the to possibility to minimise the crosstalk; ISvanda et al.l 120111 ). The 
vertical velocity is usually not evaluated from helioseismic inversion and is usually modelled 
from the horizontal components assuming the mass conservation. In our case all three flow 
components are inverted for independently (an example of the full flow map displaying 
horizontal components by arrows and vertical by colour is displayed in Fig. H]). 

To test qualitatively the sensibility of the vertical flow inversion, we utilise the 
comparison of the vertical flow with the divergence of the horizontal flow, both measured 
by time-distance helioseismology. The averaging kernels are almost identical (see Fig. [1]), 
thus the flow components are inverted for consistently. One would expect (due to the 
expected mass conservation of plasma flows and negligible horizontal variations in density) 
that divergent regions in the horizontal flow correspond to upflows in the vertical flow. Such 
behaviour is visible by eye in Fig. HJ Quantitatively the correspondence can be evaluated 
by means of correlation coefficient between Vz and d^Vx + dyVy, which is summarised in 
Tab. [3l Correlation is positive and significant. The correlation coefficient is expected to be 
higher, however one has to bear in mind that the horizontal and vertical flow components 
have different signal-to-noise ratios. That is estimated as RMS(fQ,)/crQ, and has values 4.0 
for horizontal, 1.1 for vertical components. Consequently, the signal-to-noise ratio of the 
horizontal divergence has value of 1.6. Thus one cannot expect perfect correspondence. 
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4.1. Comparison with surface Dopplergram 

An additional surface measurement of the plasma flow exists, which can also be used 
to independently verify the results of the time-distance inversion for flows. It is the direct 
measurement of the Doppler component of velocity v^, high-cadence series of which is 
used to mea s ure th e travel times. Such comparison was already done in the past e.g. by 



Gizon et al. 



( I2OOOI ). In order to compare with the Doppler component, one has to project 
the velocity vector in the Cartesian system {v:i:,Vy,Vz) to line-of-sight component fios and 
two complementary transversal components using the heliographic coordinates of 

each point of the flow map which are known from the deflnition of the Postel projection. 
The transformation equations are given in a matrix form for illustration. 



1 
COS ip{r) — simper) 
sin(y9(r) cos (p{r) 



cos'd{r) — sin'i?(r) 

1 
sin'(9(r) cos'd{r) 



Vx{r) 
Vz{r) 



(6) 



where 'i9(r) = l{r) — Iq and if{r) = b{r) — bo are the heliocentric coordinates, computed 
from the Carrington coordinates I and b. 

The reference Dopplergram for each of four datacubes was derived as a temporal 
average of all high-cadence Dopplergrams over 12 hours. As in the case of the comparison 
to the horizontal flow derived from the granule tracking discussed in Section [3], the direct 
Dopplergram is obtained with a different effective spatial resolution. Thus it must contain 
information on smaller scales than the results from the flow inversion, which acts as noise 
and naturally decreases the correlation with the projected line-of-sight velocity. Thus, as 
in Section [3l we convolved the averaged Dopplergram with the horizontal section of the 
inversion averaging kernel. An example of three comparable maps (Dopplergram, spatially 
averaged Dopplergram, and a line-of-sight velocity projection) is displayed in Fig. [5l 
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The Pearson correlation coefficient between the latter two for each of the four datacubes 
investigatec 



Gizonetal 



is giv en in Table [31 The correlation is positive and significant. Note that 
(|2000[ ) mention slightly higher correlation between the Dopplergram and a 
line-of-sight projection (correlation coefficient of 0.7). The lower correlation coefficient in 
our study originate from the region at the disc centre, where the weak (and quite noisy) 
vertical component Vz of the flow dominates the projected line-of-sight velocity. When 
the Vz is excluded (set to zero) from the line-of-sight velocity projection, the correlation 
coefficient is by 0.05 higher for all four datacubes. When the centre-of-the-disc region is 
omitted (central 70 Mm), the correlation coefficient increases by additional 0.03. 



4.2. Density scale height 

From a direct comparison between the horizontal and vertical components of the flow 
we estimated the density scale height. Let us start from separating the density p and 
velocity vector v into time-averaged components p and v and fluctuations p' and v'. Thus 
the continuity equation now reads 



^ + V-{(p + p')(v + v')} = 0. (7) 
By taking the time average of Eq. ([7]) and assuming that p' = and v' = we obtain 



V- Vlnp + V-v,, + — ^ + =V-pV = (8) 

oz p 



where v/, is the average horizontal velocity. Finally we get 



(9 In p „ _ dvz 



t^.^^ + V-v/, = -— ^-=V-pV (9) 

oz oz p 
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where we assumed that the averaged density depends only on depth. We see that if the 
right-hand side is neghgible, then the horizontal divergence and the vertical velocity are 
correlated. Moreover, the slope of the correlation gives the scale height of the density 
Hp = — [din p/dz]~^ . Since horizontal velocities and vertical ones are from the same layers, 
the correlation slope gives a weighted average of the density scale height at the sampled 
depth. 

Doing the exercise, the le ast-squares fit of the data by a linear law using the known 



error levels in both variables fiYork et al. 



20041 ) gives an estimate for the density scale 



height in the near-surface laye rs, which we find around 180 km. Such a value is of the same 



order as the 125 km found by 



November! (119941 ) from a study of the supergranulation. B y 



19961 ) 



simply taking the density scale height from Model S ( jChristensen-Dalsgaard et al. 
and integrating it over the inversion averaging kernel, we get a value of 535 km, which is 
three times larger than the scale height determined from the inversion. However, one should 
note that the density scale height determinantion is a sanity check of inversion consistency. 
Additionally, due to the presence of additional terms in Eq. ([9]), one cannot expect a better 
than an order of magnitude agreement with the solar model. 



5. Conclusion 



We measured surface horizontal flows using both coherent structure tracking and 
3-D time-distance inversions. The measurements from the two techniques are in good 
agreement. Our future goals will focus on studies of plasma dynam ics especially on studies 



of solar convection (see 



Gizon &: Birch 



2012 



Hanasoge et al. 



2012) 
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Cube No. 


Time range 


bo 


lo 


No. missing frames 


1 


12 May 2012 00:01:30-12:01:30 TAI 


-3.00 


324.94 


12 


2 


12 May 2012 12:01:30-24:01:30 TAI 


-3.00 


318.33 


6 


3 


13 May 2012 00:01:30-12:01:30 TAI 


-3.00 


311.72 


6 


4 


13 May 2012 12:01:30-24:01:30 TAI 


-3.00 


305.11 


6 



Table 1: Datacubes used in the comparative analysis: temporal coverage and Carrington 
coordinates (/q for longitude and Bq for latitude) of the central point. Also number of missing 
frames (out of 961) are given for a curious reader. 




z [Mm] 

Fig. 1. — Cuts through averaging kernel of inversion for horizontal flow {vx, identical for Vy) 
and vertical flow (f^) as a function of depth. 



V (est) 



V (est) * K 



V (time-distance) 




50 100 150 200 250 50 100 150 200 250 50 100 150 200 250 

X [lyim] X [Mm] x [Mm] 



Fig. 2. — Example maps of horizontal flow velocities [vx upper row and Vy bottom row) 
measured from datacube No. 3. Displayed are (from left) CST-tracked velocity, granule 
motion velocity convolved with the inversion averaging kernel, estimate of the flow component 
from the inversion. 



Cube No. 






slope(t;^,cst,^^a',td) 




1 


0.78 (0.60) 


0.75 (0.57) 


0.97 


0.94 


2 


0.78 (0.59) 


0.73 (0.56) 


0.94 


0.89 


3 


0.82 (0.63) 


0.76 (0.56) 


0.96 


0.93 


4 


0.79 (0.61) 


0.73 (0.56) 


0.95 


0.85 



Table 2: Statistical values for comparison of results of the two methods. We give the corre- 
lation coefficient CC of the horizontal flow components and the slope of the least-squares fit 
to the fiow estimates obtained by the two methods. 
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v^(cst) * K v^(cst) * K 



Fig. 3. — Direct comparison of values of flow estimates from the two methods for independent 
points in the field of view. The line with slope of unity represents the location where all 
points should lie in case both flow estimates are identical. 
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E, 



Fig. 4. — An example of the full-vector flow map, horizontal components are displayed by 
(reference arrow indicates flow 250 ms^^), vertical component is presented by colours. 



arrows 
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Cube No. 


CC(f^, + dyVy) 


CC(t;f,i;ios) 


1 


0.50 


0.59 


2 


0.47 


0.57 


3 


0.49 


0.64 


4 


0.52 


0.51 



Table 3: Correlation between vertical flow Vz and the horizontal divergence of the horizontal 
flow and correlation between the measured Doppler velocity and the vector velocity field 
projected into the line of sight. 
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50 100 150 200 250 
X [Mm] 



Fig. 5. — Comparison of the measured Doppler component of the velocity (left), the measured 
Dopplergram convolved with the averaging kernel (middle), and the time-distance vector flow 
estimate projected to the line-of-sight (right). One sees a striking correspondence between 
the middle and right panels. 



